function res=thetapuntod(t,X)
    global b;
    global I;
    global Ki;

    res(3)=(Ki)*(thetaR(t)-X(1));
    res(2)=(-b/I)*X(2)+(1/I)*X(3);
    res(1)= X(2);
end